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Abstract 

The validity of the renormalization group approach for large N is 
clarified by using the vector model as an example. An exact difference 
equation is obtained which relates free energies for neighboring values 
of N . The reparametrization freedom in field space provides infinitely 
many identities which reduce the infinite dimensional coupling constant 
space to that of finite dimensions. The effective beta functions give exact 
values for the fixed points and the susceptibility exponents. 
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The two-dimensional quantum gravity is important to study string theories and 
is useful as a toy model for higher dimensional quantum gravity. The matrix model 
gives a discretized version of the two-dimensional quantum gravity. Exact solutions 
of the matrix model have been obtained for two-dimensional quantum grav- 

ity coupled to conformal matter with central charge c < 1. The result can also be 
understood by means of the continuum approach ||-[§[]. Although several exact 
solutions of the matrix model have been obtained, it is worth studying approxima- 
tion schemes which enable us to calculate critical coupling constants and critical 
exponents for unsolved matrix models, especially for c > 1. In order to make use of 
such a scheme, we need to make sure that the approximation method gives correct 
results for the exactly solved cases. 

Recently Brezin and Zinn- Justin have proposed a renormalization group ap- 
proach to the matrix model 0. They drew an analogy between the fixed point of 
the renormalization group flow and the double scaling limit of the matrix model @. 
They observed that a change N — > iV + 5N can be compensated by a change of 
coupling constants g — > g + Sg in order to give the same continuum physics. They 
needed to enlarge the coupling constant space as in the Wilson's renormalization 
group approach || . Consequences of their approach have been examined by several 
groups IJ. A similar approach has been advocated previously for the 1/N expan- 



sion in a somewhat different context [nj . In the case of the one-matrix model with 
c < 1, Brezin and Zinn- Justin obtained reasonable results for the fixed point and the 
susceptibility exponents in the first nontrivial approximation. In order to demon- 
strate the validity of the renormalization group approach, however, one should show 
that the systematic improvement of their approximate evaluation converges to the 
correct result. To this end, it is important to study a model in which a renormal- 
ization group equation can be derived exactly, even if it is simpler than the matrix 
model. The vector model has been proposed for a discretized one-dimensional quan- 
tum gravity, in the same way as the matrix model for a discretized two-dimensional 
quantum gravity! 11 HH1- 



The purpose of our paper is to clarify the validity and the meaning of the renor- 
malization group approach by considering the vector model. For the vector model, 
we obtain an exact difference equation which relates the free energy — log Z^-2{g) 
to the free energy — log Z^(g — 25 g) with slightly different values of coupling con- 
stants. We find that these coupling constant shifts 5gk are of order 1/N and occur in 
infinitely many coupling constants. We also obtain infinitely many identities which 
express the freedom to reparametrize the field space. By using these identities, we 
can rewrite the flow in the infinite dimensional coupling constant space as an effec- 
tive flow in the space of finite number of coupling constants. The resulting effective 
beta function determines the fixed points and the susceptibility exponents. The 
inhomogeneous term in the effective renormalization group equation serves to fix 
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non-universal (analytic) terms of the free energy. To illustrate the procedure by an 
explicit example, we analyze in detail the cases of one and two coupling constants 
which give the first two multicritical points m = 2 and 3. We obtain the fixed points 
and the susceptibility exponents which are in complete agreement with the exact 
results. 

We first recall the renormalization group approach for the matrix model. The 
partition function Zj^{g) of the matrix model with a single coupling constant g for 
the quartic interaction is defined by an integral over aniVxiV hermitian matrix $ 



z N (g) 



-N ( -$ 2 + 



The matrix model gives the random "triangulation" of two-dimensional surfaces [|TJ. 
The 1/N expansion of the free energy F(N,g) 



F(N,g) = -—log Z N (g) = £ N-*f h (g) 

h=0 



(2) 



distinguishes the contributions fh from the surface with h handles. In the double 
scaling limit 



N 



oo, 



9^9* 



with Nn (g — g^) fixed, 



(3) 



the singular part of the free energy satisfies the scaling law 
exponent 70 + hjx linear in the number of handles h 



with the susceptibility 



fh(g) = {9 - 9*) 



2-70-71^ 



a h + 



F(N,g) = (g-g^fiN-tg-gJ). 



(4) 



The functional form of / can be determined by a nonlinear differential equation 
(string equation) p[. The continuum limit is achieved at the critical point g —>■ g* 

where the average number of triangles diverges. The double scaling limit suggests 
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that we may draw an analogy between and the momentum cut-off A . 

It has been proposed that the free energy F(N, g) of the matrix model satisfies 
the following renormalization group equation M 



AT A - f3(g)^L + 1 (g) 

dN Hyy! dg iyyi 



F(N,g)=r(g), 



(5) 
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where [3(g) is called the beta function. The anomalous dimension and the inhomo- 
geneous term are denoted as j(g) and r(g) respectively. A fixed point g* is given 
by a zero of the beta function. The exponents of the double scaling limit can be 
given by the derivative of the beta function at the fixed point 71 = 2/ ( 5'(g*), 70 = 
2 — 7(g*)//? / (fi'*)- As a practical method to obtain the beta function, Brezin and 
Zinn- Justin proposed an approximate evaluation of the beta function by integrat- 
ing over a part of degrees of freedom of the N x N matrix (7], [10|. They found 
a reasonable result at the first nontrivial order. We have computed higher orders 
of the matrix model following their approximation methods, but we find no signs 
of improvement as we proceed to higher orders. Before discussing the result of the 
matrix model, we find it more illuminating to study the case of the vector model 
where we can clarify the situation more fully. 



The partition function of the O(N) symmetric vector model is given by 



Zw(g) 



d 



N 



>exp 



fe=i 



2k 



i2\fc 



(6) 



where <fi is an N dimensional real vector [11] ||14|| . Here we introduce infinitely 
many coupling constants gj., since we need all possible induced interactions after a 
renormalization group transformation even if we start with a few coupling constants 
only. The 1/N expansion of the logarithm of the partition function gives contribu- 
tions from h loops as terms with N l ~~ h . The vector model has the double scaling 
limit iV — > 00 with N l ^ x (g — g*) fixed, where the singular part of the free energy 
satisfies the scaling law [ITT] - [|14[ 



-log 



Z N ( 9l = l,g k = 0(k>2)) 



sing 



EN 1 -" (9 

h=0 



^) 2 -^V + ---. (7) 



One should note that the power of 1/N and the definition of the susceptibility expo- 
nents are slightly different from the matrix model. Therefore the relation between 
these susceptibility exponents and the derivative of the beta function becomes 



7i 



7o 



j(g*) 



In the spirit of the approximation method of ref.[ 
iV + 1-th component a of the vector 0jy+i 



we can integrate over the 



<Pn+i 
Zn+i{9) 



(<f>N,<x), 



d ' (fi^da exp 



-(N 



(9) 
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Neglecting higher order terms in 1/N, we obtain 



(g) 



Z N+1 ( 9l = 1, g k = (k > 2)) 



J d N (f)Nexp 



gk + 5g k 



k=i 



2k 



where the shifts 5g k of the coupling constants are found to be 



k=\ 



k 



\k=l 



9_k x k-i 
9i 



k=l ^ 



+ 



(10) 



Now we shall show that the shift 5g k can be evaluated exactly in the vector model 
without using approximation methods such as described above. To demonstrate that 
the result is exact, we start with the partition function Z^_ 2 (g) . After integrating 
over angular coordinates in R^ -2 , we perform a partial integration in the radial 
coordinate x = <b 2 



(g) 



TC 2 



' N 



IT 2 



1 Wo 



dx x N ' 2 ~ 2 



exp 



■(^-2)E 



r(f- 



- dxx^-'ly: 



gkx 



.k-l 



k=l 



exp 



2k 



n 

-(^-2)E|- fc 



.(12) 



Identifying the right hand side with ((N — 2)gx/2ir)Zx(g — 25 g), we obtain an exact 
difference equation for the logarithm of the partition function 



[(-log Z N (g)) - (-logZ N _ 2 (g))] - log 



(N-2) 9l 



2tt 



-[(-log Z N (g - 25g)) - (-\ogZ N (g))\. 



(13) 



We find that the shifts 5g k of the coupling constants are exactly identical to the result 
(|TT|) of the approximate evaluation. We would like to stress that no approximation 
is employed to obtain eq. (|13|). Therefore we can infer that the vector model offers 
an example to justify the approximate evaluation method of ref.0. 

In the A" — > oo limit, we can obtain a differential equation from the exact 
difference equation (|13|) 



OO Q 

ft dg k 



log Z N (g)). 



(14) 



One can bring the quadratic term in the potential to the standard form 2 /2 since 
g>i can be absorbed by a rescaling gi<f) 2 — > 2 . We have 



Z N (gi,g2, g*, ■■■) = gi N/2 z N (i } g 2 /gf, g 3 /g'i, ■■■)■ 



(15) 
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Therefore it is convenient to use the rescaled coupling constants g k together with g\ 
as independent coupling constants 



9k = 9hl 'g\- 

We shall define the free energy for the vector model 

F{N,g) 



1 i v ( \ 1 1 N 9i 
logZATig) log . 

N 2 & 2tt 



(16) 



(17) 



If we use gi, g 2 , #3, • • • as independent coupling constants, we find from the rescaling 
identity ( |T5D that the partition function depends on gi only through the factor 
g l . Therefore, the free energy F(N,g) is independent of g\ and is a function of 
h (k > 2) only. 

We denote the partial derivatives with respect to gi,g 2 ,g s , ■ • • by \g, and those 
with g h g 2 , g 3 , ■ ■ ■ by \ g 



_d_ 

dgi 



_d_ 

dgi 



g k=2 I 



d 




d 




1 d 


dg k 


9 


dg k 


a 


9i dg k 



Thus we obtain a renormalization group equation for the free energy F 



N 



d_ 

dN 



k =2 \9i 9i ) 09k 



+ 1 



F(N,~g) = NS gi 



2gi 



Eq.(^) shows that the anomalous dimension is given by 

l(9) = 1, 

which implies a relation between two susceptibility exponents 

70 + 71 = 2. 

We read off the beta functions in the rescaled coupling constants as 

A® = N (% - . 

V 9i 9i J 

These beta functions can be evaluated explicitly using eq . (p|) as 

h(g) = -^2-3^ + 2^3, 

h{g) = -2^3 + gl ~ Shh + 3^4, 



(3 A (g) = -3<? 4 - 9 2 + 4g 2 g 3 - 2^ 3 - 8g 2 g A + Ag 



5- 



(19) 



(20) 



(21) 



(22) 



(23) 
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Let us investigate the simultaneous zero of the beta functions fa = (k > 2) in 
the spirit of ref.0. The first condition fa — gives g^ in terms of §2- The second 
condition fa = determines in terms of #2 and g^, and hence in terms of gi- As 
can be seen in eq. (p3[) , fa{g) turns out to be a sum of &<7fc+i and a polynomial in 
< Therefore we find that there always exists a solution of fa = (k > 2) for 
each given value of the coupling constant g 2 . In the following, we shall show that 
this strange result of the apparent existence of the one-parameter family of fixed 
points is due to a misinterpretation of the renormalization group flow. 

The key observation is the ambiguity to identify the renormalization group flow 
in the coupling constant space. Though the above equation (|I9|) seems to describe 
a renormalization group flow in the infinite dimensional coupling constant space, 
the direction of the flow is in fact ambiguous because all the differential operators 
(d/dg~k) are not linearly independent as we will see shortly. To see this, note that 
the partition function (^) is invariant under reparametrizations of the integration 
variable <fi. Since the model is O(N) invariant, we can obtain new informations 
only from reparametrizations of the radial coordinate x = 2 . Since the radial 
coordinate should take values on the half real line, we consider the most general 
reparametrization which keeps the integration range [0, 00) 



x 




(24) 



where e/s are infinitesimal parameters. Substituting 
with respect to e/s, we obtain a family of identities 



) in ( |T2"D and differentiating 



LjZ N (g) = 

L; 



d 



=j+l d 9t 



1+ N) 3 0g- 



(25) 
(26) 



The differential operators Lj constitute half of the Virasoro algebra. We can show 
that this algebra is identical with the one found in [ Tl\ and |TJ| . 

It is more useful to use the rescaled coupling constants g^. The reparametrization 
identity corresponding to Eq is nothing but the infinitesimal form of the rescaling 
identity ( |T5| ) and reads 

^- F(N,~g) = 0. (27) 
dgx ~ g 

In terms of the rescaled coupling constants, the reparametrization identity corre- 
sponding to £1 is given as 



iV + 2 
2N 



+ E 

1=2 



1 + a7> + M £ 4 



F(N,~g) = o, 



(2f 
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The reparametrization identities corresponding to Sj is given by 



2j\ . d ~ _ d 



F(iV,fif) = 0. 



(29) 



We see that derivatives of the free energy in terms of infinitely many coupling con- 
stants gk are related by infinitely many reparametrization identities. Thus one can 
expect that only a finite number of derivatives are linearly independent. The exact 
difference equation JIB] ) combined with the reparametrization identities (p8|) and 
( |29| ) constitute the complete set of equations to characterize the renormalization 
group flow in our approach. 

To illustrate the use of the reparametrization identities, we shall first take the 
case of a single coupling constant. Let us consider a point in the coupling constant 
space 

(01, fa, h, h, ■■■) = (9i,h, o, 0, . . .). (30) 

At this point, the j-th identity relates dF/dgj + 2 to dF/d§j + i except for the case 
of j = 1 where an extra constant term is present. Therefore we can express 
dF/dgk (k > 3) in terms of dF/dc/2 by solving these reparametrization identi- 
ties recursively in the one- dimensional subspace ( |3Q|) of coupling constants. It is 
most convenient to organize the solution in powers of 1/N. We find explicitly at 
leading order 

Bk-^ + Rk + OiN" 1 ), (31) 



dg 



k uy 2 



where 




A = 4# 2 + 1. (34) 

These solutions at leading order are sufficient to eliminate dF/dgk (k > 3) in 
favor of dF/dc/2 in (|19|). Thus we obtain a renormalization group equation with the 
effective beta function P e ^(c/2) and the inhomogeneous term r(g 2 ) 



F(N,~92)^r(g 2 ), (35) 
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where 
3 efT/~ 



k=2 

1-A) 2 t 



A 



log 



1 + VA 



3^ 2 + 4 



2^A 
32 / 



A. 



(36) 



N 1 ^ 
r (#2) = —Sgi - - + 22 R A(92, 9k = 0,k>3) 

1 1 k=2 



2VA 



log 



i - Va" 



log 



i - Va" 



- + - log - + - 2 + - H #2 + ~ + 

4 2 5 2 3 V 4/ 15 V 4/ 



(37) 



The effective beta function /3 (#2) exhibits a zero at #2 — —1/4. Furthermore, 
we can calculate the susceptibility exponents from the derivative of the effective beta 
function by using (§) and ( p0| ) 



Ti 



P'(~9 



2*1 



3 
2' 



7o 



7(^2* 



1 

2' 



(3? 



The fixed point and the susceptibility exponent are in complete agreement with the 
exact results for the m = 2 critical point of the vector model corresponding to pure 

no 

gravity pT| [O]. The beta function f3 (§2) has also a trivial fixed point at #2 — 0, 

np 

which is ultraviolet unstable since <9/3 /<9g2 ( #2 = 0) = — 1 < 0. 

We can extract the complete information from the renormalization group flow 
in our approach, namely the exact difference equation and the reparametrization 
identities by a systematic expansion of the free energy in powers of 1/N 



F(N,~92) = T, N ~ h Mh)- 

h=0 



(39) 



The free energy at leading order fo satisfies an ordinary differential equation exactly, 

(40) 



fo(~92)-P eS (~92)^(~92)=r(~g2). 
092 



where the effective beta function and the inhomogeneous term are those given in 
eqs.([36|) and (p7|). We see immediately that the general solution is given by a sum of 
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an arbitrary multiple of the solution of the homogeneous equation and a particular 
solution of the inhomogeneous equation. Since both the effective beta function 
P e ^(d2) and the inhomogeneous term r(g 2 ) are analytic in g 2 around the fixed point 
g 2 = —1/4, the singular behaviour of fo comes from the solution of the homogeneous 
equation. It is important to notice that the singular term corresponding to the 
continuum physics (the so-called universal term) is specified by the beta function 
alone. Their normalization, however, cannot be obtained from the renormalization 
group equation. The analytic contributions are determined by the inhomogeneous 
term and the effective beta function. 

To obtain the full information for the free energy up to the order N , we should 
write down the solution of the recursive reparametrization identities (^) and ( [29|) 
up to the order N~ h explicitly. By inserting the solution into the exact difference 
equation (p~3| ) and by expanding it up to the power N~ h , we obtain an ordinary 
differential equation for fh 



(1 - h)f h {g 2 ) ~ /3 eff fe)§^ 2 ) = r h (g 2 ). (41) 

og 2 

We see that the effective beta function is common to all h, whereas the inhomoge- 
neous terms depend on h. Therefore we find the singular part of the free energy 
is determined by the beta function up to a normalization 

fhifa) = falsing + fh(h) analytic 

MgiUng = (h - h*) 2 -^ lh a h + ■■■. (42) 

If we sum over the contributions from various h, we find that the renormalization 
group equation determines the combinations of variables appropriate to define the 
double scaling limit. On the other hand, the functional form of the scaled variable 
is undetermined corresponding to the undetermined normalization factor for the 
singular terms of each h 

oo oo 

F(N,~g) = E N- h f h (~g 2 )sin 9 = E N ~\~92 ~ ~92*) 2 - J °- Jlh a h 

h=0 h=0 

= (~92 - ~92*) 2 -"°f(N l ^(g 2 - g 2 *)) sing . (43) 



It has been known in the exact solution of the vector model that a nonlinear 
differential equation determines the functional form / ( iV 1 / 71 (#2 — 92*)) of the 

V / sing 

scaled variables. The equation is called the string equation, or the L_i Virasoro 
constraint ||, The above consideration means that the complete set of our 

renormalization group flow equations, namely the exact difference equation and 
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the reparametrization identities, does not give the L_i constraint for the universal 
singular terms. This is in accord with our objective of the renormalization group 
approach: to obtain at least fixed points and critical exponents even if the exact 
solution is not available. On the other hand, our equations give informations on the 
non-universal analytic terms. 

We can extend the analysis to the more general situation of finitely many coupling 
constants. Let us take 

0i = l, g k ^0 (2<k<m), g k = (k>m + l). (44) 

In this subspace, each reparametrization identity involves only finite number of 
derivatives. Moreover, the j-th identity relates derivatives in terms of coupling con- 
stants gj + i,(jj + 2, • • • , (jj+m- Therefore we can solve this identity to rewrite dFj dgj +m 
in terms of dF/dg k (j + 1 < k < j + m — 1). By successively using these identi- 
ties, we find that there are precisely the necessary number of recursion relations to 
express dF/dgk (k > m + 1) in terms of dF/dgk (2 < k < m). Consequently, we 
can reduce the renormalization group equation effectively in the space of the finite 
number of coupling constants g k (2 < k < m) 



F(N, g 2 ,..., g m ) = r(g 2 , ...,g m )- (45) 



The m-th multicritical point should be obtained as a simultaneous zero of all the 
beta functions (3^ = ■ ■ ■ = (3^ = 0. The susceptibility exponent is given by 
eigenvalues of the matrix of derivatives of beta functions f^- = dfif^/dgj at the 
fixed point, which is an (m — 1) x (m — 1) real matrix. 

To illustrate the multi-coupling case explicitly, we take the case of two coupling 
constants m = 3. We can solve the reparametrization identities recursively to leading 
order in l/N, and combine Pk(g2, 03, 9k = {k > 4)) to obtain the effective beta 
functions 

B eS ( ~ ~ ) = -2a/?7 (p 3 - 7 3 ) log a + cyclic a + (3 + 7 

12 {92,93> ( a/ 3 + /5 7 + 7a) 2 («-/3)(/?-7)(7-«) «/3 + /57 + 7«' 

R eS ( ~ ~ \ = 3a P^ (P 2 ~ 7 2 ) log a + cyclic 2 

p ^ 92 ' 93 > ( af3 + (3l + ia) 2 ( Q ,_ /3 )( / 5_ 7 )( 7 _ Q; ) af3 + (3l + ia > y > 

where a, (3, 7 are the three roots of the cubic equation 

g 3 x 3 + g 2 x 2 + x - 1 = 0. (48) 

rr rr 

We find three fixed points as the simultaneous zeros of (3% n and/5| n . All of them 
turn out to be on the m = 2 critical line 

~ 93 = -l-[2 + 9~g 2 ±2(l + 3~g 2 f 2 ]. (49) 
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The first fixed point is at (5% #3) = (— 1/3, 1/27). This result agrees with the exact 
value of the coupling constants at the m = 3 critical point. We have evaluated the 

rr 

derivative matrix Qij = d/3f /dcjj by expanding the effective beta function around 
the fixed point 

nq(g2 = -l,h = ±;) = { % IV (50) 



3 u 27 



60 



One of the eigenvalues of the derivative matrix is 3/4 which gives the exact suscep- 
tibility exponent 71 = 4/3, while the other eigenvalue is 1/2 which gives an analytic 
term. The corresponding eigenvectors also agree with the exact solution. As we 
stressed before, we can regard A 1 / 71 as an ultraviolet cut-off A 2 in taking the double 
scaling limit. Since both eigenvalues are positive, the m = 3 fixed point is ultraviolet 
stable. 

The second fixed point is at g 2 = — 1/4 and g 3 = 0. We obtain the derivative 
matrix % = dpf S /dgj at (g 2 ,g 3 ) = (-|,0) 

%(s2=-~ 23=oWjj jy , (51) 

which has eigenvalues 2/3 and —1/2. The eigenvalue 2/3 correctly gives the m = 2 
susceptibility exponent 71 = 3/2. The fixed point is attractive in the direction 
(6g 2 ,5g 3 ) oc (1,0) and repulsive in the direction (Sg 2 ,5g 3 ) oc (—2,1), which are the 
eigenvectors for the eigenvalues 2/3 and —1/2, respectively. The repulsive direction 
is tangential to the m = 2 critical line (|49 ) at the fixed point (g 2 ,g 3 ) = ( — 1/4, 0). 



The third fixed point is the trivial fixed point at the origin (#2, #3) = (0,0). 

rr 

Actually we find that the effective beta function /3f (<72>£te — 0) along the g 2 axis 

rr 

is identical to the effective beta function /3f (#2) in eq. (|36]) obtained in the single 

rr ^ 

coupling constant case. This is because /3f (#2, #3 = 0) = along the g 2 axis. 



The critical line ( |4T3| ) is a trajectory of the renormalization group flow. From 
the trivial fixed point (^2,^3) = (0,0), the flow goes to the m = 3 fixed point 
(^2,^3) = ( — l/3, 1/27) along the critical line as N — > 00. Starting from a point on 
the critical line at the left of the m = 2 nontrivial fixed point (^2^3) = ( — 1/4, 0), 
the flow also goes to the m = 3 fixed point. On the other hand, the flow goes 
to infinity if it starts from a point on the critical line at the right of the m = 2 
nontrivial fixed point. We notice that all three fixed points are isolated zeros of the 
effective beta functions, since the derivative matrix is non-degenerate at these fixed 
points. Explicit evaluation shows that the effective beta function does not vanish at 
any other point on the m = 2 critical line. 

As we described earlier, a naive application of the renormalization group ap- 
proach to matrix models has given us results which show no improvement as we 
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go to higher orders. In view of our analysis of the exact result on the flow of the 
coupling constants for the vector model, one should not be surprised by these re- 
sults on matrix models. We should look for an effective renormalization group flow 
by taking account of the reparametrization freedom in the field space. Although 
the above results are derived exactly only for the vector model, we have a similar 
ambiguity to identify the renormalization group flow in the coupling constant space 
for matrix models. For instance, there are identities which form a representation of 



Virasoro algebra in matrix models [15]. There has been some work to interpret the 



Virasoro constraints as reparametrization identities in the field space in the matrix 



model []T6|. We are trying to devise a way to implement our idea of the effective 
renormalization group flow through the application of the reparametrization iden- 
tities to matrix models in order to make the renormalization group approach more 
useful. Work along this line is in progress. 

We thank B. Durhuus and T. Hara for an illuminating discussion and M. Hirsch 
for a careful reading of the manuscript. This work is supported in part by Grant- 
in-Aid for Scientific Research (S.H.) and Grant-in-Aid for Scientific Research for 
Priority Areas (No. 04245211) (N.S.) from the Ministry of Education, Science and 
Culture. 
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